Lattice expansion from isotope substitution in the Iron based superconductors 
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We estimate the changes in lattice parameters due to iron isotope substitution in the iron-based 
high temperature superconductors using a variational calculation based on the anharmonicity of the 
Fe-As(Se) bond. For BaFe2As2 and 57 Fe to 54 Fe substitution we find a c-axis expansion of 1 • 10 _3 A 
and in-plane a-axis contraction of 1 • 10~ 5 A at 250K in good agreement with experimental values on 
(Ba,K)Fe2As2 for which a negative isotope exponent a — — < was found [P.M. Shirage et 

al, Phys. Rev. Lett. 103, 257003 (2009)]. For FeSe in contrast we find an expansion for both a and 
c lattice parameters. We discuss the relevance of these isotope induced changes of lattice parameters 
to the isotope effect in light of the well established sensitivity to the Fe-As(Se) coordination. 
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There is evidence to suggest that the iron-based 
superconductors^ are unconventional in the sense that 
superconductivity is not caused by electron-phonon 
interactions.^ Pointing in this direction is the close prox- 
imity, or even coexistence, of magnetic order or strong 
magnetic fluctuations^ as well as the results from density 
functional theory (DFT) of mundane electron-phonon 
interactions!^ 

With this in mind, the observations of significant ef- 
fects on T c with iron isotope substitution is remarkable. 
In the first study of the isotope effect by Liu et al. an iso- 
tope exponent a = — ^j" ^ ~ *-*-35 for iron substitution 
in K doped BaFe2As2 and F doped SmFeAsO was founcP 
suggestive of a more or less conventional electron-phonon 
pairing mechanism. In subsequent work by Shirage et al. 
on nominally the same (K,Ba)Fe2Fe 2 a startling negative 
isotope exponent a w —0.18 was foundPSuch a negative 
isotope effect is highly unconventional while at the same 
time the large magnitude clearly indicate a non-negligible 
influence of the lattice on superconductivity. The few 
other studies of the isotope effect on T c in these materi- 
als have ranged from tiny, a ~ —0.02, in optimally doped 
SmFeAsOi-jJ 3 to large positive, a « 0.81, in FeSei-^P 
It is of course possible given the still early stages of the 
study of these materials, that not all of these results rep- 
resent genuine isotope effects. 

At the same time, the underdoped cuprate supercon- 
ductors, which are the prototypical unconventional su- 
perconductors, also display large and doping dependent 
isotope exponents.^ Particularly intriguing is the strong 
enhancement around 1/8'th doping in La2- £C Sr :c Cu04 
where charge and spin density wave order which also 
couples to the lattice is known to influence superconduc- 
tivity. Optimally doped cuprates, with the maximum 
transition temperatures, on the other hand show little or 
no isotope effect, a possible reason why this phenomenon 
has been less studied than many other exotic phenomena 
in these materials. 

It is clearly an important but apparently largely unre- 
solved issue what is the cause of isotope effects in these 
high temperature superconductors. One clue of the rich 



interplay between lattice and electrons is the tetrago- 
nal to orthorhombic structural phase transitions and the 
accompanying el ectro nic orders that are found in both 
types of systems! 10 * 11 * In this work we consider one par- 
ticular aspect of isotope substitution in the iron based 
superconductors, namely the effect on lattice parameters 
that follows from anharmonicity in any real material. An 
important question is to compare these calculations to ex- 
perimentally measured isotope induced changes in lattice 
parameters. 

In previous worlP^l a B ls , Fe dominated phonon (en- 
ergy near 20meV) was studied by Raman spectroscopy 
and the temperature dependence of the energy and life- 
time was found to be well reproduced within a calculation 
of the phonon self energy due to phonon-phonon inter- 
actions induced by assuming a cubic anharmonic correc- 
tion to the Fe-As bond potential. A rough estimate of 
the isotope induced changes of lattice parameters were 
also made based on the expansion of an isolated Fe-As 
bond. Substitution with a lighter mass gives an expan- 
sion of the bond which was taken to imply a proportional 
expansion of both a- and c-axis lattice parameters. The 
calculation with an isolated bond may give an order of 
magnitude estimate but to get a more accurate descrip- 
tion it is important to place the anharmonic bond in the 
appropriate crystal environment of an harmonic lattice. 
This is done in the present work. The main result is that 
for realistic parameters the expansion of the Fe-As bond 
will occur predominantly in the softer c-axis, accompa- 
nied with either expansion or contraction of the in-plane 
bond lengths. For BaFe2As2 we find changes that are 
quantitatively in line with those measured in the isotope 
measurements by Shirage et al. (5 C = 3 ± 1 • 10 _3 A and 
5 a = 0±1-10~ 3 A) which gives independent support that 
those experiments were indeed done on properly isotope 
substituted samplesP^ 

To explain the calculation in a transparent fashion we 
will start by considering the one-dimensional anharmonic 
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chain. The Hamiltonian for the problem is 

H = T,^ + l fa+i - r J ~ °) 2 - § fo+i - ^ - «) 3 

with and tj the momentum and position operators of 
ion j and a the equilibrium interatomic distance for the 
harmonic problem g = 0. The aim being to calculate the 
change in equilibrium distance 8 =< rj + i — rj — a > for 
the anharmonic problem g =/= 0. The anharmonic oscilla- 
tor is classic problem which cannot be solved exactly but 
for small g it can be estimated variationally^ in terms of 

the Hamiltonian H' = £V i^T + f ( r :?'+i~ r j — ( a +^)) 2 an d 
minimizing the free energy with respect to the variational 
parameter S. Thus minimizing the thermal expectation 
value < H — W >h', using rj = j(a + S) + Xj, such 
that H' is now the standard Harmonic chain expressed 
in terms of the deviations Xj from equilibrium positions, 
and expanding to second order in g gives 

2 < H - H' >„,= kS 2 - gS < (x J+1 - x f > H > (2) 
j 

Since k > this has a minimum given by S = 
|^E ? 2 (1 - cosg) < x q x- q > where x q = 
(1/%/A) e^Xj. (With N the number of sites in the 
chain.) The expansion is thus clearly related to the mean 
square displacement of the atom, or in the zero temper- 
ature limit the zero point motion. For the infinite chain, 
solving the harmonic chain in the standard fashion, we 
find 

6= ^nk^ /^Vl-cosgcoth^ (3) 

with ujq = v2(i — cos q)k/m. The actual isotope substi- 
tution induced expansion is thus given by the difference 
in anharmonic expansions for the two masses mo and mi , 

Ai so — ^mi ^mo ' (^0 

In the zero temperature limit Eq. [3] reduces to S(T — 
0) = y=pv2, with the corresponding isotope shift 

A iso (T = 0) = ^(l/^/mi- showing that sub- 

stitution with a lighter mass (mi < mo) gives an expan- 
sion (for g > 0). Although the anharmonic contribution 
to the lattice spacing, i.e. the thermal expansion, Eq. 
[3j grows monotonically with increasing temperature, the 
actual isotope shift Ai SO (T) decreases with temperature 
from a maximum at T — 0, as the mass enters both in 
the prefactor l/y / ro and in the phonon dispersion uj q . 

For a three dimensional crystal the variational expres- 
sion, Eq. [2j generalizes to 

2{H-H') H ,= (5) 

^2 %(% • 4) 2 - 9ij(sij ■ 4) < (% • & - ^)) 2 >H> 

<ij> 



where the sum runs over all bonds i,j with harmonic 
(anharmonic) coupling constant fcy (gij), where §ij is a 
unit vector in the bond direction, and where Xj, is the 
displacement of ion i with the variational expansion 
of the bond. The expectation value is then evaluated in 
terms of the full phonon spectrum of the crystal. We 
have performed the calculation for the BaFe2As2 crystal 
assuming harmonic bonds except for the nearest neigh- 
bor Fe-As bond which is given a cubic anharmonicity. 
We assume a uniform expansion of the unit cell and con- 
sider only the tetragonal crystal structure. There are 
thus two variational parameters, the expansions 5 a and 
S c of the in-plane and out of plane lattice constants a 
and c, respectively. Our harmonic model includes a mini- 
mal number of parameters needed to reproduce a realistic 
phonon spectrum with harmonic couplings as taken from 
the fit of the harmonic (Born- von Karman) model to ex- 
perimental data of the phonon spectrum from inelastic 
X-ray scattering .151 We use the nearest neighbor Fe-As 
bond (k =4.37 eV/A 2 ), Fe-Fe (k =0.81 eV/A 2 ), and Ba- 
As (k =0.50 eV/A 2 ) and the equal height As- As (k =0.25 
eV/A 2 ). The anharmonicity is g = 103eV/A 3 for the Fe- 
As bondPln Fi gure[l]is shown the calculated changes of 
lattice parameters Ai SOya and Ai sa ^ from replacing 57 Fe 
with 54 Fe for two different temperatures. (We use 57 Fe 
instead of the naturally more abundant 56 Fe in order to 
compare to the experiments in Ref[B) where such samples 
were prepared.) As expected based on the solution to 
the ID chain, the isotope shift decreases with temper- 
ature. The results are plotted as a function of the As 
elevation angle above the Fe plane 8, see Fig. [3j keeping 
the bare lattice parameters a and c fixed. (29 + a = 180°, 
with a the Fe-As-Fe tetrahedral angle.) The results are 
consistent with the expansion of the anharmonic Fe-As 
bond, which for small angles corresponds to expansion 
in the plane and for larger (physical) angles correspond 
to expansion in the c-direction. The actual magnitudes 
of expansions and contractions are related to the larger 
stiffness of the lattice in-plane than out-of-plane. 

We have also considered FeSe with the expectation 
that the missing spacer layer would make the crystal less 
prone to expand in the c-direction. Here we were not able 
to find the relevant spring constants in the literatur e but 
made a rough fit to LDA calculated phonon spectraP^ 7 -' 
In contrast to the Bal22 crystal we did find an expansion 
in the plane of similar magnitude to the expansion in c 
as shown in Figure [2] Experimentally, a very large iso- 
tope effect a s» 0.8 was found in FeSei_z together with 
a contraction in c and expansion in a of similar mag- 
nitudes ~ 3 — 5 • 10 _4 AP Given the uncertainty in our 
parameters for FeSe we tried making the interlayer cou- 
plings stronger but found that we had to make interlayer 
Fe-Se bond stronger than the corresponding intralayer 
bond in order to get a contraction in c together with an 
expansion in a. This may indicate that there are some 
problems with the isotope substituted samples used in 
Ref. |S1 In general, we have not found contractions in c 
for any reasonably physical set of parameters. However, 
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FIG. 1: Calculated expansion of lattice parameters of 
BaFe2As2 due to isotope replacement of 57 Fe with 54 Fe as 
a function of the elevation angle of As above the Fe plane, 9. 
Expansion of c is dashed and a is solid. The bare values are 
c w 13A, a « 3.9A and 9 « 35°P 
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FIG. 2: Calculated expansion of lattice parameters of FeSe 
due to isotope replacement of 57 Fe with 54 Fe as a function of 
the elevation angle of Se above the Fe plane, 9. Expansion 
of c is dashed and a is solid. The bare values are c w 5.5A, 
a w 3.8A and 9 « 40.5° P 



possible sources of error in addition to the uncertainty in 
parameters could include only considering the undoped 
systems as well as the assumption of a uniform expansion 
of the unit cell. 

Clearly, isotope induced changes of lattice parame- 
ters are small, as estimated in this work and as seen 
experimentally 5 *, at most of the order of a few times 
10 _3 A or of the order of a 0.1% change. Based on such 




FIG. 3: Sketch of calculated changes in lattice parameters 
from substitution with a lighter Fe isotope for BaFe2As2 (left) 
and FeSe (right). 



small changes of the crystal parameters it is difficult to 
see how to generate the large isotope shifts on the su- 
perconducting transition temperature in models for non- 
phonon mediated superconductivity.^ However, it has 
been argued in the past as well as recently that small 
changes in lattice para meters can indeed reproduce large 
isotope exponent d 18 * 19 l In addition, a characteristic fea- 
ture of the iron based superconductors is that they are 
extremely sensitive to small changes in the crystal struc- 
ture such as the anion (As,Se,Te) height above the Fe 
planeP It was suggested by Khasanov et al. that the 
varying isotope exponents could be reconciled by having 
two contributions, one from the standard isotope shift 
due to Debye frequency shift and a separate contribution 
from changes in the anion height where for different sys- 
tems isotope replacements would increase or decrease the 
height W' Although an interesting suggestion, our calcula- 
tions do not support this picture, we always find a c-axis 
expansion for the lighter isotope for realistic parameters. 

We believe these calculations can provide a benchmark 
for experiments on isotope substituted samples by mea- 
suring the changes in lattice parameters. If and how these 
small changes are enough to cause large isotope shifts of 
T c or whether genuine effects of phonons on the pairing 
need to be considered remains to be explored further. 

Finally we would also point out how the present cal- 
culations complement calculations of lattice parameters 
using density functional theory. DFT calculations cannot 
directly address the effects of isotope substitution as the 
atomic mass does not enter a calculation for the ground 
state energy of the electronic charge distribution. Har- 
monic coupling constants and the corresponding phonon 
spectrum are of course readily available and needed as 
input to our calculation. Although the strength of an- 
harmonicity of the bonds would also be possible to esti- 
mate we have not found any such studies on these mate- 
rials, and have instead relied on an independent estimate 
of anharmonicity based on a many body calculations of 
the phonon self-energy and the corresponding tempera- 
ture dependence of the energy. It is however possible 
to estimate the isotope induced expansion from DFT 
by including the zero point phonon energy in the total 
energy.^ Although significantly more demanding than 
the present study such calculations would provide an in- 
teresting comparison. Another natural continuation of 
the present work on the isotope induced expansion is to 
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explore the related thermal expansiorP3 which may also 
provide an independent experimental probe of the isotope 
substituted systems. 
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